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(57) The method comprises a histogram of the radi- 
ographed image, i.e., a real histogram (step 2), a math- 
ematical model of the image chain and the object ob- 
tained by calibration. The mathematical model of the im- 
age chain and the object and a set of parameters of ac- 
quisition, of the detector, of the positioner and of the ob- 
ject are used to determine (step 1 ) two values of gray 
level, gray min and gray max, delimiting a useful zone 
and the part below gray min and the part above gray 
max are suppressed (step 3) in the real histogram, i.e., 
a limited histogram. A set of rules is applied (step 4) to 
the limited histogram in order to determine the level of 
brightness (WL) and the contrast (WW) is obtained from 
the brightness and possibly from one or more parame- 
ters selected by the user or fixed a priori. 
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Description 

[0001 ] This invention relates to a method for automat- 
ic determination of the brightness and contrast of a dig- 
itized radiographic image of an object. 
[0002] It finds a particularly interesting application in 
the medical field in which fine analyses of radiographic 
images are made, notably mammography. 
[0003] In the medical field, diagnoses generally are 
based on the study of a radiographic image. The 
progress of digital systems now makes it possible to 
vary the characteristics of the image in order to best ef- 
fect a diagnosis. Thus, it is evident that the quality of the 
image is an essential point. 

[0004] As is well known, image acquisition introduces 
numerous parameters related to the acquisition chain: 
These include the parameters of acquisition (target ma- 
terial of the anode of the X-ray tube [track], material and 
thickness of the filter at the output of the X-ray tube [fil- 
ter], voltage applied between cathode and anode of the 
X-ray tube [kV], heating current of the filament of the 
cathode of the X-ray tube [mA], exposure time, ...), the 
parameters of the positioner (incidence of the view, en- 
largement factor, type of compression pad used, thick- 
ness of the object, force of compression, ...) and the pa- 
rameters of the digital detector (relation of gray levels 
to flux of X-rays captured ...). In addition to this, the pa- 
rameters of the object, in particular its composition, are 
introduced. This set of parameters makes it possible to 
obtain an image that is in fact composed of pixels of 
shades of gray. Quantification of the image is obtained 
on a given number of gray levels. When the image is 
displayed on a screen or imprinted on a film using a giv- 
en range of gray levels, the contrast perceived may be 
considerably different among images of the same object 
acquired with different sets of parameters of acquisition 
or among images of different objects obtained with the 
same set of parameters. In addition, selection of the 
range of gray level to be considered for the screen is of 
great importance in order to give an acceptable percep- 
tion of contrast in the image. Thus, once the image has 
been obtained, the user interactively modifies the bright- 
ness (WL: window level) and the contrast (WW: window 
width) of this image in order to adjust the range of gray 
levels in it so as to obtain a correct perception of con- 
trast. 

[0005] Generally, for more than 16000 different gray 
levels, it is a difficult and lengthy task to select the cor- 
rect brightness and contrast manually. 
[0006] An embodiment of the invention permits auto- 
matically determining the value of the brightness (WL) 
and to derive the contrast (WW) from it. 
[0007] An embodiment of the invention therefore pro- 
poses a method of digital radiographic image acquisition 
of an object with automatic adjustment of the parame- 
ters of visualization (brightness, contrast). The method 
is adaptable, since the wishes of the user may be taken 
into account. 



[0008] This method comprises the obtaining of two 
gray levels, min gray and max gray, revealing a zone of * 
gray levels in which the brightness (WL) is capable of 
being determined, a set of rules for calculation of the 

5 value WL, and from this, lastly, the derivation of a range 
of gray levels for the contrast WW. 
[0009] Other advantages and characteristics of the in- 
vention will appear upon examination of the detailed de- 
scription of a mode of implementation, in no way limita- 

10 tive, and of the accompanying drawings, wherein: 

Figure 1 is a flowchart of a mode of implementation 
of the method according to an embodiment of the 
invention; 

15 

Figure 2 illustrates the limited histogram with re- 
spect to the real histogram; 

Figure 3 illustrates a way of selecting WL; 

20 

Figure 4 describes the method of obtaining an in- 
variant contrast. 

[0010] According to a general embodiment of the in- 
25 vention, a first step, segmentation, is carried out. The 
input data of this first step are the parameters of acqui- 
sition, the parameters of the digital detector, the param- 
eters of positioning and the parameters of the object. All 
of these parameters are input in a pre-established math- 
30 ematical model of the image chain and of the object. 
[0011] The model makes it possible to determine a 
gray min level and a gray max level. The part surround- 
ed by these two levels is the useful zone of gray levels 
within which the value of the level of brightness WL will 
35 be selected. 

[0012] This step corresponds to segmentation of the 
image by double thresholding so that only the points of 
the image having gray min level and gray max level are 
preserved. 

40 [001 3] A second step is carried out in parallel with the 
preceding step, it performs the real radiographic acqui- 
sition of the object. From the image obtained, a histo- 
gram is drawn up, i.e., a real histogram. 
[0014] The third step introduces gray min, gray max 

45 and the real histogram. Their correlation results in sup- 
pression in the real histogram of the part below the level 
of min gray and of the part above the level of max gray. 
The result of the third step is an interval of gray levels 
limiting the useful zone, i.e., limited histogram. The val- 

50 ue of the brightness WL within this interval then remains 
to be determined. 

[0015] Accordingly, the fourth step concerns the de- 
termination of WL, which may be done by means of a 
set of simple rules known to persons skilled in the art. 
55 This set of rules makes it possible to obtain a value WL. 
However, when the user performs numerous mammog- 
raphies, it is necessary to improve efficiency by intro- 
ducing a mathematical model of the histogram. One 
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model is described by Mr. Jean Lierrard, LSD/AAP ref- 
erence 98030 Technical Note GE Medical Systems, Au- 
gust 1 998. It considers a compressed breast, the shape 
of which is modeled, for example, with the aid of simple 
geometric forms. When data of the breast to be mam- 
mographed are introduced, a histogram model charac- 
terized by a maximum at the level of radiologic thickness 
corresponding to a breast tissue is obtained. 
[0016] Correlation of the mathematical model of the 
histogram so determined and the limited histogram 
makes it possible to determine the value of a radiologic 
thickness characterizing the breast. 
[0017] The limited histogram, thus refined, will be sub- 
ject to a set of rules permitting determination of the 
brightness WL. 

[0018] The fifth step concerns determination of the 
contrast WW. This value may be obtained in a manner 
known to those skilled in the art by introducing the value 
WL and any parameters known a priori and dependent 
upon the user. 

[0019] According to one mode of implementation of 
the invention, the contrast WW may alternatively be de- 
termined independently of WL. In this case, the coeffi- 
cient of mean attenuation (u. mean), dependent upon the 
spectrum, is used. In a general way, the WW thus ob- 
tained is a function of the spectrum, because it is pro- 
portional to u. mean. 

[0020] Obtaining WL and WW (independently of WL) 
is equivalent in fact to an auto-contrast operation that 
may be used in a method of perception of contrasts in 
invariant thickness. 

[0021] An embodiment of this method makes the per- 
ception of contrasts independent of the conditions of ac- 
quisition and of the object. For two objects of different 
thicknesses, it is desired that regardless of the condi- 
tions of acquisition, the contrast perceived on the image 
should remain faithful to their difference in real thick- 
nesses. 

[0022] In other words, it is desired that a given thick- 
ness should always represent the same perception of 
contrast, regardless of what the object and the condi- 
tions of acquisition are. 

[0023] In effect, in the course of image acquisition, the 
physical elements cause the spectrum to convert the re- 
al thickness of the object into radiologic thickness, thus 
resulting in a perception of unfaithful contrast. This is an 
effect of exponential attenuation. The effect is compen- 
sated for by introducing into an image chain a change- 
of-space step in order to annul the exponential attenu- 
ation due to interaction of the X-rays with the object by 
employing a modified logarithmic function. The change 
of space makes it possible to leave the exponential 
space to pass into the space of radiologic thicknesses. 
[0024] Likewise introduced is a visualization step 
making it possible to pass from the space of radiologic 
thicknesses to the space of real thicknesses. This 
change of space is possible because WW is proportion- 
al to the coefficient of mean linear attenuation of the ob- 



ject studied and independent of WL. 
[0025] Although the invention is not limited, the meth- 
od is applicable to the automatic determination of the 
brightness WL and of the contrast WW for a mammog- 
5 raphy. 

[0026] As shown in Figure 1 , the first step employs a 
mathematical model of the image chain and the object 
with, in input data, the following parameters: 

10 - thickness of the compressed breast and parame- 
ters of the positioner (incidence of the view, enlarge- 
ment factor, type of compression pad used, thick- 
ness of the object, force of compression, ...) as pa- 
rameters of the positioner, 

15 

parameters of the detector (relation between the 
flux of X-rays received on the detector and the gray 
levels of the image produced, ...) 

20 - parameters of acquisition (track, filter, kV, mAs, ...) 

parameters of the object (mechanical thickness of 
the breast, minimum u. min and maximum u. max val- 
ues of the coefficient of linear attenuation of the 
25 object, ...). 

[0027] Since the breast is composed principally of fi- 
brous and adipose tissues, if there is no information on 
the composition of the breast, u. min and u. max may be 

30 estimated by making two extreme assumptions. 

[0028] u, min corresponds to the coefficient of linear 
attenuation of the least attenuant tissues of the object 
(adipose tissues for the breast) for the energies of the 
X-ray spectrum determined by the parameters of acqui- 

35 sition. 

[0029] u. max may be estimated in two ways: on the 
one hand like jx min, by considering that u. max corre- 
sponds to the coefficient of linear attenuation of the most 
absorbent tissues of the object (fibers for the breast) for 

40 the energies of the X-ray spectrum determined by the 
parameters of acquisition; on the other hand, in a more 
precise manner, from a mathematical model of the im- 
age chain, the mechanical thickness of the compressed 
breast, the parameters of acquisition and a quantity of 

45 photons obtained following preexposure performed on 
a zone of maximum density (which makes it possible to 
estimate the value of the coefficient of linear attenuation 
corresponding to the most attenuant zone of the object). 
[0030] The set of parameters introduced into the 

50 mathematical model of the image chain and the object 
makes it possible to have two values of gray levels, gray 
min and gray max, at the output (Fig. 2). These two val- 
ues, as a matter of fact, delimit the useful zone, which 
is a zone of gray levels relating really to the breast. In 

55 effect, owing to the two extreme values of the composi- 
tion of the breast, u, min and u, max, a zone delimited by 
two extreme values, gray min and gray max, has been 
obtained, outside of which the gray levels do not corre- 
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spond to the breast. More precisely, the part of gray lev- 
els below gray min corresponds to objects more atten- 
uant than the object of interest, and the part of gray lev- 
els above gray max corresponds to the bottom of the 
image. This step is a segmentation step because it s 
makes it possible to delimit the useful zone. 
[0031] Then a correlation is effected between the two 
values gray min and gray max and a histogram obtained 
from the radiographic image of the breast, i.e., a real 
histogram. More precisely, the part below gray min and 10 
the part above gray max are eliminated so as to pre- 
serve only the useful zone: limited histogram (Fig. 2). 
[0032] The brightness WL is a value included in the 
useful zone and may be obtained in a variety of ways. 
[0033] One manner of obtaining WL is the application is 
of a set of pre-established rules to the limited histogram. 
A set of rules may include: 

determination of the gray level corresponding to the 
maximum of the limrted histogram 20 

preservation of a quantity of x% (typically 95%) of 
occurrences of the limited histogram on the right of 
the maximum, and likewise x% of occurrences on 
the left of the maximum: a reconstructed histogram 2s 
is thus obtained 

determination of WL as median value of the recon- 
structed histogram (Fig. 3). 

30 

[0034] It is alternatively possible to obtain WL with 
better precision by preceding the set-of-rules step by a 
correlation step. This step introduces a mathematical 
model of the histogram in which: 

35 

the shape of the breast is a cylinder generated by 
rotation about an axis of a rectangle of which one 
of the short sides is closed by a semicircle equal in 
diameter to the length of this short side 

40 

the composition of the breast is homogeneous, for 
example 100% fat 

a h istogram is established that corresponds to prob- 
ability as a function of the radiologic thickness of the *s 
breast 

the maximum of the histogram obtained represents 
the maximum thickness of the breast which, multi- 
plied by the coefficient of attenuation, gives the so 
maximum radiologic thickness corresponding to the 
adipose tissue (fat). 

[0035] Correlation of the two histograms (mathemat- 
ical model and limited histogram) makes it possible to 55 
determine in the limrted histogram the value of the max- 
imum radiologic thickness corresponding to the adipose 
tissue in the breast. 
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[0036] This then makes it possible to determine the 
values of radiologic thickness of the various compo- 
nents of the breast. 

[0037] This correlation is performed by employing a 
method of minimization of errors between two functions 
such as, for example, the method of least squares. 
[0038] An appropriate set of rules can then be applied 
in order to determine WL. For example, WL = ccE, with 
E representing the value of the radiologic thickness ob- 
tained from the mathematical model of the histogram. 
[0039] The following step corresponds to the determi- 
nation of WW by using WL Thus, WW is obtained from 
a function introducing WL 

and possibly other parameters, in particular a pa- 
rameter G (Fig. 3), which is selected by the user 
[0040] This parameter G therefore makes the method 
adaptable to each user. 

WW = g(WL, G) 

■ g being a function which, from WL and G, deter- 
mines WW, which thus represents a range of gray levels 
about WL. 

[0041] It is alternatively possible to determine WW in- 
dependently of WL. 

[0042] First, u. mean is determined from information 
drawn from the limited histogram. As an example, the u. 
corresponding to the median value of the limited histo- 
gram may be taken as \x mean. 
[0043] Then, a law introducing a constant Cte is used 
in order to derive WW from it: 

WW - Cte u. (spectrum) 

[0044] This relation is true in a mono-energetic case 
but in a general way, WW is a function of the spectrum 

WW = f (spectrum) 

[0045] The operation of auto-contrast has thus been 
performed, since WL and WW have been established. 
[0046] The way in which perception of the contrast of 
a difference in thickness remains invariant regardless of 
the means of acquisition and the objects will now be ex- 
plained with reference to Figure 4. 
[0047] The imposition of X-rays on the object results 
in an exponential atten uation of the intensity I at the level 
of the image: 

I = lo exp (-J u.dl) 

[0048] lo is a constant, 1 represents an infinitesimal 
magnitude that corresponds to a distance along the path 
connecting the focus of the X-ray and the detector. 
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[0049] J u.dl represents the radiologic thickness for a 
given zone of an object. 

[0050] It is this quantity that interests us. To obtain it, 
a pre-LUT (look-up table) operation is performed, mak- 
ing it possible to offset the exponential attenuation by 
using a modified logarithmic function. It is called modi- 
fied because the lowest gray levels are converted ac- 
cording to a linear rotation when the logarithmic function 
is progressively introduced for the other gray levels. * 
[0051] We thus find ourselves in the space of radio- 
logic thicknesses in which a radiologic thickness may 
be given the notation nH, with H the real thickness. 
[0052] Lastly, a change of space must be made to get 
back into the space of real thicknesses. 
[0053] Auto-contrast furnishes us with the values WW 
and WL, which are introduced at the level of a visuali- 
zation lut. The said visualization lut makes it possible to 
eliminate u. In the case of a mono-energetic image, this 
operation amounts to division by u~ 
[0054] This operation is possible because WW is pro- 
portional tO [L 

[0055] The result thus obtained may be introduced in- 
to the DICOM standard visualization system (Grayscale 
Standard Display Function, Supplement 28), familiar to 
persons skilled in the art, in order to visualize the image 
(Fig- 4). 

Claims 

1 . A method of acquisition of a digital radiographic im- 
age of an object, having a histogram of the radio- 
graphed image (real histogram) and a mathematical 
model of the image chain and the object obtained 
by calibration, comprising the steps of: 

a) providing the mathematical model of the im- 
age chain and the object and of a set of param- 
eters of acquisition, of the detector, of the posi- 
tioner and of the object, to determine two values 
of gray level, gray min and gray max, delimiting 
a useful zone 

b) suppressing the part below gray min and the 
part above gray max in the histogram (limited 
histogram); 

c) applying a set of rules to the limited histo- 
gram in order to determine the level of bright- 
ness WL; and 

d) obtaining the contrast WW from WL and pos- 
sibly from one or more parameters selected by 
the user or fixed a priori. 

2. The method according to claim 1 , comprising the 
step of estimating a coefficient of minimum linear 
attenuation (u. min) of the object from known values 



of the coefficient of linear attenuation of the least 
attenuant tissues of the object (adipose tissues for 
the breast) for the energies of the X-ray spectrum 
determined by the parameters of acquisition and 
5 makes it possible, with the parameters of the detec- 
tor, the parameters of the positioner, the parameters 
of acquisition and the parameters of the object, 
through a mathematical model of the image chain 
and of the object, to determine gray min. 

10 

3. The method according to claim 1 , comprising the 
step of estimating a coefficient of maximum linear 
attenuation (u. max) of the object either from known 
values of the coefficient of linear attenuation of the 

is most attenuant tissues of the object (fibrous tissue 
for the breast) for the energies of the X-ray spec- 
trum determined by the parameters of acquisition 
or, in a more precise manner, from the mechanical 
thickness of the compressed breast, the parame- 

20 ters of acquisition, the mathematical model of the 
image chain and the object and from a quantity of 
photons obtained following a preexposure per- 
formed on a zone of maximum density of the object, 
and makes it possible, with the parameters of the 

25 detector, the parameters of the positioner, the pa- 
rameters of acquisition and the parameters of the 
object, to determine gray min through a mathemat- 
ical model of the image chain and the object. 

30 4. The method according to any one of the preceding 
claims, wherein the limitation of the histogram intro- 
duces two gray levels (gray min and gray max) sur- 
rounding a zone of gray levels within which the 
brightness WL is determined. 

35 

5. The method according to any one of the preceding 
claims, wherein a correlation is made between the 
limited histogram and a mathematical model of the 
histogram in order to determine the value of a radi- 

40 ologic thickness characterizing the object. 

6. The method according to claim 5, wherein the cor- 
relation between the limited histogram and the 
mathematical model of the histogram is effected by 

45 applying a method of minimization of errors be- 
tween two functions. 

7. A method of acquisition of a digital radiographic im- 
age of an object, having an acquisition chain com- 

50 prises a step of compensation for the effect of ex- 
ponential attenuation of the radiation by using a 
modified logarithmic function, so that perception of 
the contrast (WW) of a given difference of thickness 
remains invariant regardless of the means of acqui- 

55 sition. 

8. The method according to claim 7, wherein the 
brightness (WL) and (WW) obtained from a coeffi- 
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cient of mean attenuation of the image independ- 
ently of the brightness, are introduced in a visuali- 
zation step so that the signals coming from the com- 
pensation step (space of radiologic thicknesses) 
are sized in the space of real thicknesses. s 
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